Dynamics at T ^ in the half-integer isotropic high spin molecules 
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We investigate the dynamical spin-spin auto-correlation function of the isotropic high spin molecules 
CrCue (S = 9/2), CrNie (S = 15/2) and CrMne (S = 27/2 ), using magnetization, ^SR and NMR 
measurements. We find that the field autocorrelation time r of the molecule's spin at zero and low 
fields is nearly temperature independent as T ^ 50 mK. The high temperatures r is very different 
between the molecules. Surprisingly, it is identical (~ 10 nsec) at low temperature. This suggests 
that r is governed by hyperfine interactions. 
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High spin molecules (HSM) consist of clusters of metal 
ions, they are ordered in a crystal lattice, and coupled 
by Heisenberg ferromagnetic or antiferromagnetic inter- 
actions with coupling constant J, only between spins Si 
in the molecule. At low temperatures fc^T < J these 
spins lie parallel or anti-parallel to each other, and the 
molecule is in its ground spin state, where S = Si is 
high, with quantum number S, and 25-1-1 degeneracy. At 
even lower temperatures, fcsT ^ J the degeneracy can 
be removed by additional anisotropic interactions such as 



etc. 



the uniaxial term DS^, or rhombic term E[S'^ 
Experiments on the two most famous high spin molecules 
Mni2 [0 and Fes show that at low T the main in- 
teraction is, indeed, the uniaxial anisotropy, where up 
and down spin states Sz = ±5 = ±10 are degenerate, 
and tunneling is induced between these states by an ad- 
ditional term in the Hamiltonian that does not commute 
with Sz. This quantum tunneling of the magnetization 
(QTM) phenomenon has received considerable attention 
in recent years, focusing mainly on the additional terms 
in the Hamiltonian which are responsible for the tun- 
nehng |^|-||. Nevertheless, the theoretical picture is far 
from being clear and existing models are controversial, 
and often contradict each other ||^. In some cases, even 
qualitative understanding of the observed experimental 
data is absent Q . 

In this paper we investigate three simple HSM systems, 
which are isotropic {D, E = 0). In these systems no tun- 
neling is observed due to the absence of the uniaxial term 
DSz- However, spin dynamics is observed even at very 
low temperatures (T = 50 mK). Therefore, the additional 
terms in the spin Hamiltonian could be probed directly 
in these systems. Such a study can highlight the role of 
phonons [|6|j^], dipolar interactions, and nuclear fluctua- 
tions 1^ in a simple setup. In addition, it could serve as 
a large S limit for isotropic models usually applied only 
for the S = 1/2 case or it could serve as a I? ^ test 
case for anisotropic HSM models ^ . 

We present an experimental investigation using three 
types of measurements: magnetization, spin lattice re- 



laxation {T^ ^) of muon spin (/^SR), and of proton 
spin (NMR). Our molecules are [Cr{(CN)Cu(tren)}e 
](C104)2i 0, [Cr{(CN)Ni(tetren)}6i(C104 )9 
and [Cr{(CN)Mn(tetren)}6] (0104)9 which we la- 

bel as CrCug, CrNig and OrMng respectively. In these 
molecules a Or (HI) ion is surrounded by six cyanide ions, 
each bonded to a Cu(II), Ni(II) or Mn(II) ion. The co- 
ordination sphere of Or and Cu/Ni/Mn can be described 
as a slightly distorted octahedral. The objective in this 
work is to find the spin-spin correlation time t{T) in the 
three systems, and to compare them. Our main findings 
are: (I) r is nearly T independent at low temperatures, 
and (II) at very low temperature r does not depend on 
5 or J in this isotropic case. 
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FIG. 1. The magnetization of CrCue, CrNie and CrMne 
as a function of the external applied field at T = 2 K, re- 
spectively. The solid lines are the S = 9/2, S — 15/2 and 
S — 27/2 Brillouin functions (see text). 

In Fig. we show the magnetization AI per molecule 
in units of /is as a function of applied field H for the 
three molecules. The data are taken at T = 2 K. In all 
cases M increases as the applied field is increased. The 
magnetization reaches a saturation value of | gi^B, ^ 
gUB and ^ gfiB CrCug, CrNig and CrMng respec- 
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tively. These saturation values are consistent with 6 Cu 
(5 = 1/2), Ni {S = 1) or Mn {S = 5/2) ions ferromag- 
netically (or anti-ferromagnetically in the case of CrMng) 
coupled only to a Cr ion of spin 3/2. The solid lines are 
Brillouin functions of the respective S. In Refs 10 -T^] 
the susceptibility was fitted to the one expected from 
the Heisenberg Hamiltonian and the values JcrCue = 77 
K, JcrNig = 24 K and JciMne = — 11 K where found. 
Therefore, the highest spin value of each molecule is well 
seperated from other spin states (a few tens of degrees K). 
High field ESR measurements (on CrNig) |l^ and sus- 
ceptibility measurements (on CrCue, CrNig and CrMng) 
found no evidence for anisotropy, namely, D ^ 0. This is 
consistent with the octahedral character of the molecules. 

In our fiSK experiments we measure the polarization 
P(t, H) of a positive muon spin implanted in the sam- 
ple, as a function of time t and magnetic field H, where 
P{0,H) = 1. The field is applied in the direction of the 
initial muon polarization. The positive muon decays to 
a positron which is emitted in the direction of the muon 
spin, and the polarization as a function of time is recon- 
structed by the detection of the emitted positrons. 

The measurements in all molecules are done at tem- 
peratures ranging from 25 mK up to 300 K, and in fields 
ranging between zero and 20 kG. These experiments were 
performed at both ISIS and PSI, exploiting the long time 
window in the first facility for slow relaxation (high T), 
and the high time resolution in the second facility for fast 
relaxation (low T). 



H=0-100GT[K]= 
■ 2 T 27 

o 5 42 
* 17 + 200 




FIG. 2. The spin polarization as a function of time, (a) At 
zero field and difi'erent temperatures, (b) At field H — 20 kG 
and different temperatures. The solid lines are fits of the data 
to square root exponential functions. 

In Fig. H(a) and (b) we present the muon spin polariza- 
tion as a function of time and for different temperatures 
in CrNig in zero field, and in _ff = 20 kG, respectively. In 
zero field, the relaxation rate increases, as the tempera- 
ture is decreased, and saturates at 5 K. The increase 
at high temperatures is caused by thermally activated 
tarnsitions between excited spin states. However, at low 
temperatures, only the ground spin state is populated. 



and only transitions within the degenerate ground state 
are possible. In contrast, at = 20 kG, and tempera- 
tures lower than ~ 17 K, the relaxation rate decreases as 
the temperature is decreased, and does not saturate. 

In Ref. we demonstrated that the magnetic field 
experienced by the muon in all molecules is dynamically 
fluctuating even at T = 50 niK. We therefore analyze our 
data using spin lattice relaxation theory. In this theory, 
the polarization of a local probe (muon or nucleus), in 
the fast fluctuation limit, is given by 



P{H, t) = {Po - Poo) exp [-t/Ti] + Po< 



(1) 



where Pq is the initial polarization, P^x 
polarization |ll6|] , 

Ti{H) = A + BH"^, 
1 
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and 7 is the probe gyromagnetic ratio. The correlation 
time T and mean square of the transverse field distribu- 
tion at the probe site in frequency units are defined 
by 



7^ (Bi(i)B±(0)) =A2 exp (-t/r) 



(3) 



The fast fluctuation limit is obeyed when rA ^ 1. 

In /iSR Poo = 0, and the muon could occupy many 
different sites in the sample, because the molecules are 
fairly large (^ 15 A diameter) and are embedded in an 
organic surrounding. As a result one must average over 
A. Using the distribution |nq 



P(A) 
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and allowing for a constant background (Pg) due to 
muons stopping outside the sample one obtains 



P(t)=Poexp -Ji/T^^ 



Ba 



(5) 



where l/r|' is the muon spin relaxation rate. This form 
is in agreement with the experimental results (see below) . 
In addition, Eq. ^ still holds, while in Eq. |^ and |2^, A 
is replaced by A*. 

The solid lines in Fig. || are fits of the data to Eq. ^ 
where Pq is a global parameter. The parameter Bg is free 
within 10% of its mean value since the high fields affect 
the positron trajectory in a manner that is reflected in 
Bg. The fit is satisfactory in all cases apart from the 
highest H and lowest T in CrNig. The relaxation rate 
l/Tl* in the different compounds, obtained from the fits 
is presented in Fig. ^ as a function of temperature for 
different values of H . As pointed out above, at low fields 
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l/T^ increases with decreasing temperatures, and satu- 
rates at low temperatures. In addition the value of 
increases as the spin of the compound is higher, as ex- 
pected from Eqs. |^ and ^, and the saturation temperature 
increases as the coupling constant J increases. This is in 
strong contrast to Mni2 and Feg where in zero 
field 1/Ti increases continuously upon cooling until the 
correlation time r becomes so long that the molecule ap- 
pears static in the muon dynamical window, and no Ti 
saturation is observed at low temperatures. 



etc.? In order to address this question we performed 
proton NMR Ti measurements. We find that in all ap- 
plied fields smaller than 20 kG the proton Ti is shorter 
than the experimental window around the peak in l/T-f'. 
Only in a field of ^ 20 kG were we able to perform the 
experiment at all temperatures. We measure Ti using a 
saturation— i — tt/2 — n pulse sequence. The proton po- 
larization recovery follows Eq. |l| with Pq = from which 
we obtained T/^. 
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FIG. 3. l/Ti as a function of temperatures for different 
external fields, measured by /xSR and by NMR (after scaUng) 
in (a) CrCue (b) CrNie and (b) CrMng. 

In Fig. U we plot the average relaxation time (iJ) 
at T ^ as a function of for all compounds, for 
fields up to 2 kG (note the axis break). We find that 
Ti^ obeys Eq. ^ This implies that the muon spin re- 
laxation is indeed due to dynamical field fluctuations, 
and that at low T the field autocorrelation can be de- 
scribed by a single correlation time as long as the ap- 
plied field is not too strong. At high fields (> 2 kG) 
we find deviations (not shown) from the linear relation 
between Ti and H^. The deviation might be due to 
the impact of the field on the spin dynamics, i.e. the 
correlation function given by Eq. ^. From the linear 
fits in_Fig. HL and taking (A*)^ = jfj,{AB)~^/^ from 



Eq. |2b| and g^, where 7,, = 85.162 MHz/kG, we find 
that for CrCug = 4.9 ± 0.9 MHz (57 ±10 G), for 
CrNie AS = 26 ± 2 MHz (305 ± 25 G), and for CrMng 
A^ = 38±2 MHz (446 ±24 G); the subscript stands for 
T — > 0. Using T = {B/Aj'^y/'^ from the same equation 
we find To = 7 ± 1 nsec for CrCue, tq = 10 ± 1 usee for 
CrNie and tq = 11 ± 1 nsec for CrMng. These values 
of Aq and tq are self consistent with the fast fluctua- 
tion limit. Most striking is the fact that all tq values are 
nearly equal. 

Although our data support a picture where the muon 
spin relaxes due to dynamically fluctuating magnetic 
flelds, they leave open the interpretation of these fluctu- 
ations. Are they due to the fluctuations of the molecular 
spins or a result of muon diffusion, muonium formation. 




FIG. 4. The saturation relaxation time as a function of 
for CrCug , CrNig and CrMne. The solid lines are linear fits 
of the data. 

Since the proton gyromagnetic ratio — 42.57 
MHz/kG is very different from that of the muon, we scale 
the NMR results and plot them together with the /iSR 
resuhs in Fig. |, for CrNig (b) and CrMug (c). The 
scaling factor C used in Fig. || is 0.6 for CrNig and 8.8 
for CrMng. It was chosen so that the scaled NMR re- 
laxation rates will agree with the /iSR rates at high T. 
After this scaling, we flnd a good agreement between the 
/iSR and NMR data at all temperatures. In flelds lower 
than 20 kG, where the NMR Ti was measurable only at 
high T, we obtained the same agreement between the two 
techniques. 

The scaling factor provides information on the ratio of 
the fleld experienced by a muon and by a proton. At high 
temperatures where Ti shows no fleld dependence one can 
assume that A ^ BH^ in Eq. 0, and therefore 1/Ti = 
A^T. Using the deflnition of A given in Eq. we can 

write C = hHll^f (T/^/rf) = (Bi) ^ (Bi)^ where 
(B^) is the RMS of the transverse fleld at the proton 
site, and ^B^)^ is the RMS of the fleld at the muon site 
in its general sense given by Eq. ^ The proximity of 
C to 1 especially in the CrNig is very encouraging and 
suggests that the muon sites are close to a proton in this 
system. Thus, we prove that in both techniques we are 
measuring the probe's spin lattice relaxation time due to 
the molecular spin fluctuations. 

Finally, we would like to obtain the correlation time 
at all temperatures. This could be calculated from 
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combined with magnetization measurements, at zero (or 
very low) fields. In zero order approximation we assume 
that (B^) is proportional to (S^) which is different from 
S{S + 1) since at temperatures fcsT ~ J, states other 
than the ground state S can be populated. Therefore in 
zero field 
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Taking 



{S') {T) 



S{S + 1) 



3kBTx{T) 



(6) 



(7) 



where N is the number of molecules, /is is Bohr magne- 
ton, kg is the Boltzman factor, and g — 2, we find 



r(T) 



S(S + l) 



3kB T,{T)Tx{T)A, 



(8) 



In the insets of Fig. ^ we present (5*^) (T), obtained 
from _ff — > DC-susceptibility measurements and Eq. |^, 
as a function of temperature for CrCug, CrNig and 
CrMng. In Fig. || we present r(T) as calculated using 
Eq. ^ for the different compounds. The T dependence of 
the correlation time t{T), unlike the muon spin lattice re- 
laxation rate, reflects the dynamics of the molecular spin 
without the T dependence of the field at the muon site. 
At T ^ 100 K there is more than an order of magnitude 
difference in r between the different molecules. As the 
temperature is lowered the correlation time in all com- 
pounds increases as the temperature is decreased, but 
reaches a common saturation value of ~ 10 nsec (within 
experimental error) at T 10 K. At this temperature 
only the ground state S is populated. In other words, 
when the HSM are formed they all have the same corre- 
lation time at low T. 



At very low temperatures one can interpret the corre- 
lation time T in terms of a broadening of the spin lev- 
els, due to interactions between a molecular spins and 
other molecular spins, or the enviroment. This inter- 
action is considered as a purturbation, the strength of 
which should be of the order of hr^^. Therefore, our re- 
sults have two major indications. First, the broadening 
in the T — > limit is not due to phonons since these die 
out exponentially with temperature at 0.05 — IK. Sec- 
ond, the broadening cannot be explained by interactions 
which are quadratic in S or have higher S dependence. 
This rules out the dipolar interaction between neighbor- 
ing molecules since in the three compounds the nearest 
neighbor distance is 15 A. Similarly, crystal field terms 
which are allowed by the octahedral symmetry (S"* or 
higher jl^) are unlikely. 

The only mechanisem suggested to date for level broad- 
ening of HSM, which is weakly (linear) S dependent is 
the hyperfine interaction between nuclear and electronic 
spins. This mechanism can account for the finite spin lat- 
tice relaxation rate at very low temperatures pC| . The 
hyperfine interactions in unisotropic high spin molecules 
were studied recently p2 23 1, and their effect on QTM is 
becoming clearer p^ , p5[ . We believe that this interaction 
also governs the fluctuations of the isotropic molecules at 
very low temperatures (T < 3 K). However, at high tem- 
perature (T > 10 K) the fluctuations are governed by 
spin-phonon interactions ||2l| ]. 
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